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Abstract 

We present a systematic analysis on the basis of ab initio calculations and many-body pertur¬ 
bation theory for clarifying the dominant fluctuation in complex charge-spin-orbital coupled 
systems. For a tight-binding multiband model obtained from the maximally-localized Wan- 
nier function analysis of the band structure by the local density approximation, we take into 
account electron correlations at the level of random phase approximation. To identify the dom¬ 
inant fluctuation, we carry out the eigenmode analysis of the generalized susceptibility that 
includes all the multiple degrees of freedom: charge, spin, and orbital. We apply this method 
to the paramagnetic metallic phase of a pyrochlore oxide Cd20s207, which shows a metal- 
insulator transition accompanied by a peculiar noncoplanar antiferromagnetic order of all-in 
all-out type. We find that the corresponding spin fluctuation is dominantly enhanced by the 
on-site Coulomb repulsions in the presence of strong spin-orbit coupling and trigonal crystal 
field splitting. Our results indicate that the combined method offers an effective tool for the 
systematic analysis of potential instabilities in strongly correlated electron materials. 

Keywords: ab initio calculation, random phase approximation, strong electron correlation, Cd 20 s 207 , 
spin fluctuation, all-in all-out 


1 Introduction 

The interplay of charge, spin, and orbital degrees of freedom under the electron correlation has 
been one of the central issues in condensed matter physics HIM- It brings about fascinating 
phenomena in a wide range of physical properties, including magnetism, optics, and transport. 
Among them, the properties induced by symmetry breaking and associated fluctuations have 
been of primary important, since they lead to a dramatic change of the equilibrium states 
and a gigantic response to an external field. Despite a long history of research, a systematic 
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theoretical study of complicated entanglement between the multiple degrees of freedom remains 
as a challenging issue in condensed matter physics. 

For clarifying the underlying mechanism in these complicated phenomena, it is helpful to 
identify which degree of freedom of electrons plays a role. An efficient way is to analyze the 
fluctuations. For a continuous phase transition, when approaching the critical point from the 
high-symmetry phase, a particular fluctuation diverges corresponding to the symmetry break¬ 
ing, which provides the information of the relevant degree of freedom in the phase transition. 
Even in the case of a first-order phase transition, the analysis of enhanced fluctuations might 
give useful information. Thus, identifying the dominant fluctuation in the high-symmetry phase 
is significantly important for the understanding of phase transition phenomena. 

In this paper, we present a systematic study of the dominant fluctuation in the charge, 
spin, and orbital coupled systems. Our method is twofold: the construction of the tight- 
binding model by ab initio calculations and the perturbation theory for electron correlations. 
In the first step, we perform band structure calculations on the basis of the relativistic local 
density approximation (LDA), and estimate the tight-binding parameters by the maximally- 
localized Wannier function (MLWF) analysis. In the second step, we take into account the 
effect of electron correlations at the level of random phase approximation (RPA), and calculate 
the generalized susceptibility that includes all charge, spin, and orbital degrees of freedom. 
To identify the dominant fluctuation, we carry out the eigenmode analysis of the generalized 
susceptibility while changing the model parameters as well as temperature. 

We apply the technique to a pyrochlore oxide Cd20s2C>7 and discuss the results in compari¬ 
son with those in the previous experimental and theoretical studies. Here, we briefly summarize 
the physical properties of this compound. The lattice structure of Cd20s2C>7 is schematically 
shown in Fig. QJa). Cd 2+ cations are nonmagnetic, while Os 5+ cations have three electrons in 
t 2 g orbitals per site, and comprise the pyrochlore lattice, as shown in Fig. [l|b). This material 
exhibits a metal-insulator transition at T c = 225 K mm- The phase transition is of second 
order, and no structural changes are observed. The transition is accompanied by a noncopla- 
nar antiferromagnetic order of the so-called all-in all-out type, as shown in Fig. mb) 0. This 
magnetic order occurs with the wave number q = 0, which does not break the symmetry of 
the lattice. On the theoretical side, the electronic structure of this material was discussed by 
ab initio calculations including relativistic effects [30. The peculiar all-in all-out type mag¬ 
netic order as well as the metal-insulator transition was explained by including the electron 
correlations by the relativistic LDA+t/ method [5j. 


2 Model and Method 

2.1 Model construction 

We start with ab initio calculations based on the density functional theory with LDA for the 
target material mm-- Cd 2 0s 2 0 7 in the current case. In the LDA calculation, we use a 
fully-relativistic ab initio computational code, QMAS (Quantum Materials Simulator) |12l . 
The relativistic effects are fully taken into account. We employ the projector augmented 
wave method m- The computational condition is the same as that in the previous study 
for Cd 2 0s 2 0 7 0. 

Based on the LDA band structure, we construct the multiband tight-binding model for the 
relevant orbitals: in the current case, the f2g orbitals of 5 d electrons in Os cations (see Fig. [2j. 
The multiband Hamiltonian is given by H = Ho + "Hint, where "Ho and H- ln t denote the one- 
body and two-body part, respectively. In order to estimate the tight-binding parameters in 
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Figure 1: (color online), (a) Schematic picture of the lattice structure of CCI 2 OS 2 O 7 . The blue 
and gray spheres denote Cd and Os, respectively, (b) The pyrochlore lattice structure composed 
of Os cations. The allows represent the magnetic pattern of the all-in all-out order below T c . The 
spins point in the local [ 111 ] directions, which connect the centers of neighboring OS 4 tetrahedra. 
All-in and all-out tetrahedra appear alternatively, resulting in the q = 0 noncoplanar magnetic 
order. 


the one-body part "Ho, we adopt the MLWF analysis of the LDA results mm, as follows. 
Here, Ho consists of three parts, namely, electron transfer Ht, crystal-field splitting Hcf , and 
spin-orbit interaction Hsoi . Ht is given by 

Rt = ^Cp.C'p' W C lrTpk C C<?P'k> (1) 

feElBZ C'pQ'p' (t 


where c^ ap k and fc denote the Fourier transform of the annihilation and creation operators, 
respectively, for the obtained MLWF V'CffR+p( r ); C, o, R , and p denote the orbital, spin, unit 
cell, and sublattice in the unit cell, respectively. The summation of the wave number k is taken 
in the first Brillouin zone for the Bravais lattice defined by {!?}. The hopping matrix element 
t(p,Cp' (k) is the Fourier transform of the overlap integral given by 

k,C K-R + P) ~ ( R ' + P’)] = j drtp* (aR+p (r) ip c > a > R > +p >{r). (2) 

For Hcf , we include the dominant trigonal distortion of each OsOg octahedra [5], as 


/° 1 1 \ 

Hcf = A x 1 0 1 , 

V 1 0/ 

where the basis is taken as {d yz ,d zx ,d xy ). Meanwhile, Hsoi 

A / 0 ~ i<Jz 
Hsoi = XI ■ s = - x \ ia z 0 


is given by 


l<7y 

0 


(3) 


(4) 


where the basis is taken as in Eq. ([3| including the spin degree of freedom explicitly, and a x , 
a y , and a z denote the Pauli matrices for spin indices. We estimate the coupling constants A 
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and A by using the MLWFs as 

A = (d xy t \n LDA \d yz t>, (5) 

\ = (dxy t |H LDA I d yz i ), (6) 

where 'H LBA is the Hamiltonian in the relativistic Kohn-Sham equation. 

For the two-body part of the Hamiltonian, we consider the on-site electron correlation 
defined by 

'Hint = C lo- 1 r C M<T 2 r c /x , o- 2 rC^a 1 r- (7) 

CpC'A 0102 r 

For simplicity, we assume that the MLWFs possess the same symmetry as the atomic orbitals. 
Then, the coupling constant is simplified into the form of 

EWv = (U — 2 Jn)5^'5yyi + </h(^oAv + (8) 

where U and Jh denote the intra-orbital Coulomb interaction and the Hund’s-rule coupling, 
respectively. Here, we assume the rotational symmetry of the Coulomb interaction; U = U' + 
2Jh and Jh = Jpair, where U’ and J pa j r are the inter-orbital Coulomb interaction and the pair 
hopping between different orbitals, respectively. 


2.2 Random phase approximation and eigenmode analysis 

We study charge, spin, and orbital fluctuations in the obtained multiband model by calculating 
the generalized susceptibility. First, we calculate the bare susceptibility from the one-body 
part of the Hamiltonian TLo, which is given by 

X l oX,a’^) = - T E E^'P',ap( fc )^pW( fc + ^)> ( Q ) 

1BZ 0Jk 

where a denotes the orbital and spin indices as a = (£«,cr a ), and T denotes the temperature. 
Here, G^} p > ap (q) is the noninteracting Green function, where q = ( q,to q ); ui q is the Matsubara 
frequencies. 

Next, we take into account the effect of two-body part T-Lmt at the level of RPA. For this 
purpose, we define the vertex function as 


^a/3p,/3'a'p' — (^CoC/3 6p>G a i d<j a a p , ^0fi0 a i + ^CaC/3'»C a 'C/3 d<T a 0 a id<Tpi0p)d pp ', 
and calculate the generalized susceptibility in the matrix form of 


X 


RPA / 


(«)= I-X {0 \q)U 


( 10 ) 


( 11 ) 


where I denotes the unit matrix. 

We can find the dominant fluctuation by studying the eigenvalues and eigenvectors of the 
matrix x RPA (( 7 ), as follows. According to the fluctuation-dissipation theorem, the generalized 
susceptibility satisfies the relation: 

d{ci pq cp pq )= E ^$ P , a 'p' P '(q)h a 'p' P '(q), (12) 

ot'fi'ft' 
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where S (A) represents the deviation of the thermal average of A induced by a generalized 
external field h a 'p' p '(q). If the external field is parallel to the /cth eigenvector of x RPA , e ap P i 
Eq. (fT^l) reads 

( clpqCppq) oc x K (q)e*p p {q), (13) 

where S K (A) represents S (A) in the Kth eigenmode, and x K (q) is the /cth eigenvalue. This 
indicates that S K (ca Pq cp pq ) also becomes parallel to . Thus, the eigenvector and eigenvalue 
correspond to the normal mode and amplitude of the fluctuation, respectively. Therefore, the 
analysis of the eigenmodes provides the information of charge, spin, and orbital components of 
the fluctuations. For example, the spin fluctuation Ssp(q) in the /cth eigenmode at the sublattice 
p is written as 


6s;(q) = S' 



where cr denotes the vector of Pauli matrices; <j = ( cr x ,a y ,a z ). 


(14) 


3 Application to CCI 2 OS 2 O 7 

We apply the method in the previous section to a pyrochlore oxide Cd 20 s 207 . The band 
structure for this compound was calculated in the previous study [9] for the experimental 
lattice structure at 180 K [5j. Note that the system does not show a structural change at the 
metal-insulator transition at 225 K. Figure Oja) shows the relativistic LDA band structure. 
The band structure near the Fermi level consists mainly of t 2 g orbitals of the Os cations, which 
are well separated from other bands. We carry out the MLWF analysis for the t 2 g bands and 
estimate the tight-binding parameters in Ho- As shown in Fig. [2ja) , the tight-binding band 
structure calculated by Ho well reproduces the LDA result. In this 5 d electron system, the 
spin-orbit interaction is relevant: A = 332 meV, which is larger than the largest matrix element 
in the nearest neighbor hopping, |t| = 178 meV. The trigonal crystal field splitting is also 
substantial: A = 96.6 meV. Figure [21(b) shows the density of states obtained by the MLWFs. 
There is a sharp peak of DOS at the Fermi level. In the previous relativistic LDA + U study, it 
was discussed that the Coulomb repulsion U suppresses the DOS and induces the all-in all-out 
type magnetic order [S]. In the following, we examine the electron fluctuation on the basis of 
the LDA band structure by using the eigenmode analysis of the generalized susceptibility. 

First, we compute the bare susceptibility x^ in Eq. J9J. Figure [3] shows all the eigenvalues 
of the static component, X^Hq) = X^(9> w = 0)- I n the calculations, we take 32 3 q points 
in the first Brillouin zone and set the temperature to 10 K. Although the temperature is 
much lower than T c = 225 K, we confirmed that the following results are qualitatively similar 
when the system is in the paramagnetic state in the vicinity of instability toward the magnetic 
ordering discussed below. The number of the eigenmodes is 144: we show the sublattice diagonal 
components of x^Hq) with three orbital and two spin components, namely, 144 = 4 x (3 x 2) 2 . 
Note that the eigenvalues at q = 0 are not available precisely in the present scheme due to the 
technical reason. As shown in Fig. [3l the largest eigenmode has a broad peak around the L 
point, while it is rather suppressed around the T point. However, as described below, this wave 
number dependence is changed considerably by electron correlations. 

Next, we compute the generalized susceptibility at the level of RPA by Eq. (fllT) . FigureQJa) 
shows the maximum eigenvalues of X RPA (?)- We show the data while changing the intra-orbital 
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Figure 2: (color online), (a) Electronic band structure for Cd 20 s 207 near the Fermi level. The 
curves denote the results obtained by the relativistic LDA calculations, while the crosses show 
those by the tight-binding Hamiltonian 'Ho obtained by the MLWF analysis, (b) The density 
of states calculated for Ho- In (a) and (b), the Fermi level is set to zero. The inset in (b) 
represents the first Brillouin zone for the pyrochlore lattice structure. 



Figure 3: Eigenvalues of the bare susceptibility calculated by using the tight-binding 

parameters obtained from the MLWFs. 
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Coulomb interaction U from 700 meV to 780 meV with keeping its ratio to the Hund’s-rule 
coupling as Jn/U = 0.2. Other parameters are taken to be the same as in Fig. [3j As shown in 
Fig. Oja), the eigenvalues around the T point are strongly enhanced as U increases. 



Figure 4: (color online), (a) U dependence of the maximum eigenvalues of X RFA {q)- We 
change the intra-orbital Coulomb interaction U, while keeping Ju/U = 0.2 and T = 10 K. (b) 
Spin fluctuations corresponding to the maximum eigenmode near the T point, at q = qo = 
(7r/16,7 t/16,7t/ 16). The calculations are done for U = 780 meV and Jh = 156 meV. The 
vertical axis represents the amplitudes of spin fluctuations in an arbitrary unit. Green, blue, 
and red bars represent the amplitudes of the x , y , and z spin components, respectively. The 
numbers 1-4 in the horizontal axis denote the sublattices, (c) Schematic visualization of the 
spin fluctuations in the maximum eigenmode in (b). 


We investigate the nature of dominant fluctuation by using the eigenmode analysis intro¬ 
duced in Sec. Here, we consider the eigenmode at qo = (7t/16, 7t/16, 7t/ 16), which is the 
smallest wave number along the T-L line. We find that, in the dominant eigenmode, spin 
fluctuations defined by Eq. (11411 have much larger amplitudes compared to charge and orbital 
components. Figure Efb) shows the x, y, and z components of the spin fluctuations at qo. We 
find that all the three components have almost similar amplitudes, indicating that the direc¬ 
tions of the spin fluctuations are roughly in the (1,1,1), (1, —1, —1), (—1,1, —1), and (—1, —1,1) 
directions at the sublattices 1, 2, 3, and 4, respectively. The spin fluctuation is of all-in all-out 
type, whose schematic visualization is shown in Fig. [2c). Thus, this peculiar spin fluctuation, 
which breaks the rotational symmetry in spin space, is caused by the electron correlation in 
the presence of the strong spin-orbit coupling and the trigonal crystalline electric field. The 
importance of the spin-orbit coupling and trigonal crystalline electric field was also pointed out 
in the previous LDA+t/ theory [9]. 


4 Discussions and Concluding Remarks 

Let us discuss our results in comparison with the experiments and the previous theory for 
Cd 2 0s 2 0 7 . In our results, the fluctuation is dominantly enhanced around the T point, as 
shown in Fig. [2a). In addition, the dominant fluctuation is in the spin channel and of all-in 
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all-out type, as shown in Figs. 0|b) andUJc). The results are consistent with the experiments 
in CCI2OS2O7: a q = 0 magnetic order with all-in all-out type was observed associated with 
the metal-insulator transition [Bj. We note that the peculiar magnetic order was reproduced 
by the relativistic LDA+C/ method, in which the electron correlation is treated at the mean- 
field level and fluctuations are neglected [9]. In this sense, our RPA study on fluctuations is 
complementary to the previous theory. 

Our method is not specific to this compound: it provides a versatile tool to study dominant 
fluctuations in a wide range of materials. We can also argue, for clarifying the underlying 
physics, how the fluctuations are affected by controlling the model parameters, not only the 
electron correlations but also the one-body tight-binding parameters. Moreover, the method 
can be more sophisticated by adopting other techniques to treat the electron correlations, such 
as the self-consistent fluctuation theory. These extensions of the model and method are left for 
future study. 

We acknowledge fruitful discussions with T. Misawa, M. Udagawa, and Y. Yamaji. A.U. 
is supported by Japan Society for the Promotion of Science through Program for Leading 
Graduate Schools (MERIT). The crystal structures are visualized by using VESTA 3 [16]. 
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